Spin edge states: an exact solution and oscillations of the spin current 
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We study the spin edge states, induced by the combined effect of spin-orbit interaction (SOI) 
and hard-wall confining potential, in a two-dimensional electron system, exposed to a perpendicular 
magnetic field. We find an exact solution of the problem and show that the spin resolved edge 
states are separated in space. The SOI generated rearrangement of the spectrum results in a 
peaked behavior of the net spin current versus the Fermi energy. The predicted oscillations of the 
spin current with a period, determined by the SOI renormalized cyclotron energy, can serve as an 
effective tool for controlling the spin motion in spintronic devices. 
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I. INTRODUCTION 

The principal importance of SOI is in its ability to link 
the electron charge and spin degrees of freedom, which 
is fertile for novel physical phenomena^ 2 -! 3 -. Unlike the 
charge, the electron spin is double-valued and identifies 
two system components, which can be separated as in the 
spin Hall effect-^ or mixed via the spin Coulomb drag££ 
There are different mechanisms, realizing SOI 1 , and the 
interplay between them produces another rich arena for 
study and potential applications in spintronics^^. 

In two dimensional electron systems (2DES) of the 
quantum Hall effect geometry, the extended edge state 
play a central role in understanding of transport phe- 
nomen a 10 ' 11 ! 12 . The suppression of backscattering and 
inter-edge-state relaxatio n 13 ' 14 i 15 make possible non- 
dissipative transport through edge channels. In the 
presence of SOI the nonlocal transport through spin- 
polarized edge channels holds promise of providing even 
richer phenomenology and greater power of electronic ap- 
plications. Another strong motivation for investigations 
of spin edge states is related to the recent experimental 
realization of the Mach-Zehnder— . 

Recently several theoretical papers have addressed 
the effect of SOI on the edge states in restricted 
2DES±Z»i2>±2 1 22i2I and along magnetic interfaces 2 ^ 3 .. All 
these works, however, find unlikely an exact analytical 
solution of the edge state problem and adopt different nu- 
merical approache o 18 ' 19 i 20 i 22 i 23 , use a parabolic confining 
potential 2 ^, which has no flat domain, or give an analyti- 
cal approximation in the limit of strong magnetic fields 1 - 
where the effective SOI coupling is small. 

Here we present an analytical solution to the spin edge 
states, induced by the combined effect of SOI and hard- 
wall confining potential. We derive an exact formula 
for the electron energy dispersions and calculate their 
spectral and transport properties. We find that due to 
SOI the spin edge states are resolved not only in the en- 
ergy but are also separated in space: an effect, which is 
not captured by the approximate approach, adopted in 
Ref. [l?]. From the energy spectrum we calculate the elec- 



tron group velocity and the average spin components. We 
find that the magnitude of spin components are not equal 
in the up and down states. In the quasibulk states the 
electron spins are mainly aligned along the magnetic field 
while near the hard-wall the spins of edge states become 
mainly perpendicular to their propagation direction. Us- 
ing these ingredients we calculate the components of the 
net spin current as a function of the Fermi energy and 
show that the SOI induced splitting results in the peaked 
behavior of the spin current. We argue that the predicted 
oscillations of the spin current with a period, determined 
by the renormalized cyclotron frequency, can serve as a 
new tool for manipulating spin currents in a controllable 
manner. The developed approach here is equally appli- 
cable to the spin edge states along magnetic interfaces in 
nonhomogeneous magnetic fields. 



II. THEORETICAL CONCEPT 

We assume that the 2DES resides in a quantum well, 
formed in the (001) plane of a zincblende semiconduc- 
tor heterostructure, and is exposed to a perpendicular 
homogeneous magnetic field, B = Bqz. The motion of 
electrons in the 2DES is confined by an infinite hard- 
wall potential, V(x) — oo for x < 0. Such a system is 
described by a two-dimensional Hamiltonian of the form 

H = H a + Hsoi + V(x) (1) 

where the Hamiltonian of free particle in a magnetic field 
is Ho = (7? 2 /2m*) f and the Rashba SOI Hamiltonian 
Hsoi — ctR (K x <j y — Tr y a x )2^, m* denotes the electron 
effective mass, ff = p — (e/c)A the kinetic momentum 
with p = —iKV. We choose the Landau gauge so that 
the components of vector potential are A = (0,xBq,0). 
The unity matrix f and the Pauli spin matrices a act in 
the spin space. It is assumed that electrons are confined 
to the lowest energy subband in the z-direction. 

Using the ansatz ^ (x, y) = e lkyV Xk y (x), we can reduce 
the two-dimensional Schrodinger equation H^ = E^> 
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to the one-dimensional problem where E is the elec- 
tron total energy and k y the electron momentum in 
y— direction. Then the electron wave function Xk y {%) in 
x— direction should satisfy the following equation 




(2) 



Xky 0) = 



where the effective potential V c g(x) = (x — X (k y )) /4 
in x— direction depends on the wave vector k y . In Eq. [2] 
we express the energy E — > (v + 1 / 2) huj b in units of 
the cyclotron energy, Hub = HeBo/m*c, and the length 
x — > xIb/V% in the magnetic length, Is = y/Hc/eBo. We 
introduce also the dimensionless SOI coupling constant 
7 = \/2aji/vB with the cyclotron velocity vb = h 2 /m*lB 
and the dimensionless coordinate of the center of orbital 
rotation X(k y ) — ^/2k y ls- Taking into account explicitly 
that the eigenstates of Eq. [5] are spinors, 
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Xlky (X ' 
X2ky {X ■ 



X{k y )) 
X{k y )) 
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we write the Schrodinger equation in the following com- 
pact form 



h + \ ( Xiky (£) 



K J \X2ky (£) 



. 



Here we introduce the following operators: 



h v 
h± 



de 
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The system of equations, obtained from Eq. 2] has to be 
solved under the boundary conditions \k v (x) — * when 
x — > and +oo. In the absence of SOI, h± — 0, the so- 
lution is given in terms of the parabolic cylindrical func- 
tions, D v (x). In the presence of SOI we search the bulk 
solution of matrix equation (U|) as Xik y (£) = ciD^ (£) and 
Xiky (£) = b-D/j-i (£) where a and b are the x independent 
spinor coefficients and fi is an arbitrary index, different 
from v. Making use of the following recurrent properties 
of the parabolic cylindrical functions 



KDp (0 = (u - m) £>„ (0 , 



h±D ll (0 



(0 



we obtain 
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where c± = b±/a±. Thus, the two independent bulk 
solutions of Eq. [4] are given by the spinor wave functions 



Xt(0=a± 



(11) 



The normalization of the wave functions 
/ d£,xly(0Xk y (0 = 1 gives the amplitude of eigen- 
states 



a± 



de(|^±(or + |cii a |^ ± -i(o| 



1-1/2 



(12) 

It is easy to see that in the limit of vanishing SOI, 7^0, 
we have a + — * and 6_ — * and recover the usual edge 
states, which are doubly degenerated with respect to the 
spin: 



xUx)~D v (x) 



aIld Xky ^ D " 



(13) 



On the other hand the solution (|11| for sufficiently 
large values of X (k y ) describes quasibulk Landau states 
so that the index fi± (z/, 7) differs only exponentially from 
the Landau index / = 0, 1, 2 . . . and the parabolic cylin- 
dric functions are given by their asymptotics Hermite 
polynomials, D t (£) = 2" 1 ' 2 exp (-£ 2 /4) Fj(Vf/2). In 
the limit of X (k y ) — > 00 taking n± (y, 7) = I in Eqs. © 
and (jlOp , one can exactly reproduce the spectrum and the 
wave functions of the bulk dispersionless Landau levels, 
renormalized by the SOI for I = 1, 2 . . ^^^^^ 



£f(7)=n±W~+Z 7 2 ^ 



1 




(14) 



(15) 



As usual the I = Landau level remains not perturbed 
by the spin-orbit coupling. 

In order to obtain the spectrum of the spin edge states 
we require vanishing of the electron wave function (jllj) at 
x = 0. The energy, E, obtained from the vanishing con- 
ditions for both spinor components at x = 0, should be 
the same. As seen, however, the different spinor compo- 
nents of the bulk solution (fTl"]) are given by the parabolic 
cylindric functions with different indices, which makes 
impossible their vanishing simultaneously. To satisfy the 
boundary conditions, we construct a linear combination 
of the two independent bulk solutions as 



A y (0 = a X U0 + PxL(0 



(16) 



and choose the coefficients a and j3 so that the compo- 
nents of the new spinor wave function ipk (£) vanish at 
x = 0. The eigenvalue problem for a and (3 has a so- 
lution if the respective determinant vanishes at x = 0. 
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FIG. 1: (Color online) The energy spectrum of spin edge states as a function (a) of the momentum k y for 7 = 0.3 and (b) of the 
effective SOI coupling 7 (or that is the same as a function of _B~ 1//2 ) for X(k y ) = 3. The dashed and solid curves correspond 
to the up and down spins, (c) The particle average position as a function of k y for 7 = 0.3. The probability density in the 
respective symbol positions is shown in the insets. 



This leads to the following exact dispersion equation 

C_IV (-X(ky)) iV_ X (-X(ky)) 

= c + D^_ (-X(k v )) (-X{k v )) (17) 

for the spin edge states with the wave functions 



D, + (0 - rIV (?) 



(18) 



Here r = D ll+ (—X(k y )) /D ll _ (—X(k y )) and a is ob- 
tained from the normalization of the wave functions. Re- 
call that the dependence on energy E = (y + 1/2) Ttiob 
manifests itself via the functions fi± (1^,7) and c± (v,~f), 
given by Eqs. © and (flQ|) . The dispersion relation (TT7|) 
is quadratic with respect to the parabolic cylindric func- 
tions, therefore for a given band index, n, it has two 
solutions, E sn (k y ), where s =| and J, corresponds to the 
spin edge states with the two spinor wave functions 



E=E, 



(19) 



III. 



ENERGY SPECTRUM OF SPIN EDGE 
STATES AND SPIN CURRENT 



Further we carry out the actual calculations of the 
spectrum of spin edge states, the average components 
of spins, and the spin current components, carried by 
the skipping orbits along the edges of 2DES. In the pres- 
ence of a perpendicular magnetic field, the efficiency of 
SOI is determined by the dimensionlcss coupling con- 
stant 7, which is inversely proportional to \J Bq. There- 
fore, the SOI effect is significant in weak magnetic fields 
unless temperature fluctuations smear magnetic quanti- 
zation effects. In such fields the Zeeman effect is small 21 
and we do not consider it here. We carry out the ac- 
tual calculations for Bq corresponding to the cyclotron 
splitting about 5 K. In GaAs with the electron effec- 
tive mass m* = 0.067mo such a cyclotron splitting is 
achieved for Bq w 0.25 T and taking the Rashba coupling 



a,R ~ 4.72 meV A, we have 7 = 0.03. For such a coupling 
the spin-orbit effects in the edge state spectrum should be 
hardly visible. The situation is favorable in InAs where 
the Rashba coupling is larger, an ~ 112.49 meV Ai De- 
spite the smaller effective mass, m* = 0.026mo, in weak 
fields about B = 0.1 T, we have 7 = 0.45. As we see 
below such a coupling results in essential modifications 
in the spectrum and transport of spin edge state, mea- 
surable in experiment. 

In Fig. la we plot the energy spectrum of spin edge 
states, E sn (k y ), as a function of momentum k y , which we 
obtain by solving the dispersion equations (|17p . It is seen 
that for a given quantum number n there are two spin re- 
solved magnetic edge states, E\_ n (k y ) and E-\ n (k y ). Both 
branches show monotonic behavior in the whole range of 
k y variation. For large positive values of k y , the energy of 
spin edge states is given by the spin-split quasibulk Lan- 
dau levels, renormalizcd by SOI, while at negative values 
of k y the spectrum describes the current carrying edge 
channels. The spin splitting of edge states increases with 
the main quantum number n. At the stronger effective 
SOI coupling the spectrum shows well pronounced anti- 
crossings. The development of the anticrossings can be 
traced clearly in Fig. lb where we calculate the energy 
of spin edge states versus 7 or, what is the same, ver- 
sus 1/\/Bq. These anticrossings in the energy spectrum 
result in additional structures of the spin current versus 
the Fermi energy. In Fig. lc we plot the average trans- 
verse position of the spin edge states from the boundary 
of 2DES as a function of their center of orbital motion, 
defined as 



(*y) 



dxx \ipk y (x - X(k y ))\ 



(20) 



It is seen that except for large positive values of k y , the 
position of skipping orbits takes spin resolved values so 
that the up and down spin edge states are separated in 
space. Notice this effect of the SOI induced spatial sepa- 
ration is missed in the approximation, adopted in Ref.ll7l. 
The differences of the probability density for different 
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FIG. 2: (Color online) (a) The x and z— components of average spins in units of h as a function of X(k v ) for the first two 
bands n = 1, 2 and for 7 = 0.3. (b) The a;— component and (c) the z— component of the net spin current (solid curve) versus 
the Fermi energy. The dashed and dotted curves plot the separate contributions to the spin current, made by the spin up and 
down states together for each band n. Insets show the absolute value of the separate up and down spin current contributions. 



spins and wave vectors of the first two bands are clearly 
shown in the insets. As seen from the lower inset, the 
probability density for different bands and spins differs 
even at large positive values of k y . In this limit, however, 
irrespective the quantum number n and the spin orien- 
tation, the particle average thickness is the same and 
varies linearly with its center of orbital motion. This is 
because in the quasibulk Landau states far from the in- 
terface, electrons oscillate symmetrically with respect to 
the guiding center, X(k y ), independent of the spin and 
band index. 

From the obtained spectrum we calculate the corre- 
sponding group velocities along y— direction, v sn (k y ) = 
dE sn (k y ) / dhky (v x = along x— direction) as well as the 
average spin components along x, z— directions 



QX.Z 

^ sn 



(ky) = 77 



dxipl (x)a x . z ip ky (x) 



E=E sn (k y ) 



(21) 

Because the transverse wave functions are real, the 
y— component of the spin vanishes identically, S^ n (k y ) = 
0. In Fig. 2a we plot (k y ) as a function of X(k y ) for 
the first two bands. At large positive values of k y when 
electrons are far from the hard- wall, the spins are mainly 
aligned along z— axes. This is because in the quasibulk 
Landau states electrons have no preferential direction in 
the (x, y)— plane of their cyclotron rotation. In the op- 
posite limit of negative k y , the edge channels are formed 
and the spins are mainly aligned in a;— direction, perpen- 
dicular to the y— direction of electron propagation. No- 
tice that due to the spin splitting the absolute values of 
the average spin components do not equal in the up and 
down states and this asymmetry becomes stronger with 
the band index n. 

In Fig 2b and 2c we calculate the x and z— components 
of the net spin current along y— direction, defined as 



JX,Z 



(22) 



We use this definition of the spin current, which is widely 
accepted and intuitive from the physical point of view. 



Notice that according to some recent suggestions 29 ' 30 , 
the standard definition of spin current (which is the one 
we use here) is a proper definition and there is no need 
for other definitions 31 i 32 ' 33 . It is seen that J y exhibits 
regular peaks as a function of the Fermi energy. Such an 
oscillatory behavior is due to the SOI induced splitting 
in the edge state spectrum. For a given quantum num- 
ber n the splitting always creates a narrow energy region 
near the peak position (see the inset in Fig lb) where 
only the | spin states contribute positively to the spin 
current. With an increase of the Fermi energy the | spin 
edge state starts to contribute negatively at the position 
of I spin Landau level. In this region the exponential 
increase of the velocity 3 ^ of the | spin edge state with 
the energy results in a sharp peak of the spin current. 
With a further increase of the Fermi energy the J, and 
| spin edge states of the next n + 1 band start to con- 
tribute similarly but with stronger amplitudes because of 
the spin-splitting enhancement with n. Thus the peaks of 
the spin current are imposed against the monotonic back- 
ground and have a period, determined by the cyclotron 
energy. The latter is renormalized in the presence of SOI, 
as seen in Fig. la in the limit of k y — > 00. 

As seen in Fig. 2c the z— spin current J y changes its 
sign, in addition to its peaked behavior: due to an inter- 
play between the average spin S* n (k y ) and the velocity 
v S n (ky), the spin current is negative near the Landau 
levels and positive between them. In this case, at high 
energies corresponding to large negative values of k y in 
each band n (cf. Fig. la), the large values of v sn (k y ) are 
compensated by the small values of S' n (k y ) (cf. Fig. 2a). 
Therefore the overall monotonic increase of J y with Ep 
becomes less pronounced. 



IV. SUMMARY 

In conclusion, we present an exact analytical solution 
to the spin edge states, induced by the combined effect of 
the Rashba SOI and of the hard- wall confining potential. 
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The exact solution of the problem allows its deeper intu- 
itive understanding and can be a strong input in studying 
the spin transport through edge channels. We find that 
due to SOI the spin edge states are resolved not only in 
the energy but are also separated in space. In the bulk 
of sample the electron spin is mainly aligned along the 
magnetic field while near the hard-wall the spin becomes 
mainly perpendicular to the edge state propagation di- 
rection. The magnitude of spin components is asymmet- 
ric in the up and down states. We show that the spin 
current components exhibit oscillations versus the Fermi 
energy. The predicted oscillations, with a period deter- 



mined by the renormalized cyclotron energy, can serve as 
an effective tool to control the spin motion in spintronic 
devices. 
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